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A local void modifies the sky distribution pattern of the cosmic microwave background 
(CMB) polarisation by gravitational lensing and produces B-modes from E-modes for an off- 
center observer. In order to see whether this effect can be used to observationally test the 
validity of the local void model, we calculate this lensing effect by solving the propagation of 
CMB polarisation along null geodesies close to the central light cone in the general Lemaitre- 
Tolman-Bondi (LTB) model perturbatively. In particular, we give general formulas for the 
correlations of E and B observed by an off-center observer and show that i5™ and are 
correlated for the same value of £, i.e., {E^B^ ) oc 5^,^/, while {E^E^ ) cc 5\ii-iii. This 
feature can be used to distinguish the gravitational lensing effect by a local void from those 
by normal shear field of galaxies. 

§1. Introduction 

Type la supernova (SNIa) observations imply an acceleration of the cosmic ex- 
pansion if general relativity is valid on cosmological scales and if the universe is 
homogeneous and isotropic on scales larger than 200Mpc. If we abandon one of 
these assumptions, however, other explanations become possible. One approach of 
such a nature is the so-called modified gravity that abandons general relativity. The 
other is the local void model, which was first proposed by Tomita^^^-* Goodwin 
et al.,^^ and Celerier,^^ independently. This model abandon the second assump- 
tion, which is often called the Cosmological Principle or the Copernican Principle, 
and assumes that we are around the center of a low density spherically symmetric 
void and that the spacetime is well described by the Lemaitre-Tolman-Bondi (LTB) 
model^)^-'.''-' In this model, the cosmic expansion rate decreases outward on each 
constant time slice, which produces an apparent acceleration of the universe when 
observed along the past light cone. Although this model violates the Cosmological 
Principle and requires an accidental situation concerning our location in the universe, 
it does not require any dark energy or a modification of gravity theory. Further, as 
far as the redshift-luminosity distance relation obtained by the SNIa observations is 
concerned, this model can reproduce the observational results with any accuracy be- 
cause it contains at least one arbitrary function of the radius (see, e.g., 8)). Actually, 
it has passed all observational tests so far. Therefore, it is of crucial importance to 
find observational tests that enable us to discriminate this void model from spatially 
homogeneous models employing dark energy or modified gravity theory, in order to 
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establish the necessity of dark energy or a modification of gravity. 

There have been proposed various tests for that purpose so far:^)'^°) CMB 
anisotropics on large angular scales, radial BAO,^^-* spectral distortions of CMB,^^^ 
the kinematic SZ effect/^) and estimates of the Hubble ratc.^^^ Among these, the 
oldest and simplest one is to observe the effect of the inhomogeneity on CMB temper- 
ature. This effect was first estimated by Alnes and Amarzguioui^^^ for a special class 
of void models. Recently, this analysis was extended to a wider class of models by 
Kodama, Saito and Ishibashi with the helps of analytic formulas for the dipole and 
quadrupole moments of the CMB temperature anisotropy for an off-center observer 
in a general spherically symmetric universe. Although this type of test provides 
a strong constraint on the allowed range of our distance from the void center, it 
is not so decisive because lower moments including the dipole and quadrupole are 
significantly affected by the cosmic variance. 

One possible way to circumvent this weakness is to extend the analysis of the 
off-center CMB anisotropics to polarisation. As is well known, in the spatially ho- 
mogeneous cosmology, inhomogeneities of the matter (galaxy) distribution modify 
the sky pattern of CMB polarisation by gravitational lensing and produce B-modes 
from E-modes^^^^^^^^^^^^.^^^ In the local void model, the central observer detects no 
such effect because of the spherical symmetry, in spite of the strong inhomogeneity 
of the model. However, an off-center observer can detect the gravitational shear field 
through the observation of B-modes. Therefore, in this paper, we calculate the grav- 
itational lensing effect on the CMB temperature and polarisation anisotropies for an 
off-center observer in the local void model and estimate the correlations among the 
temperature, the E-mode and the B-mode anisotropies. 

The paper is organised as follows. First, in the next section, we review the 
basic matters on CMB polarisation that are relevant to the present paper. Then, in 
Section 3, we perturbatively solve the geodesic equations in the LTB model to find 
the change in the propagation direction due to the gravitational lensing effect for 
null rays close to the central past light cone. This determines the shift vector on the 
sky for an off-center observer that represents the difference between the observed 
direction of a light ray and the corresponding direction of the last scattering point 
in the fiducial spatially homogeneous model in a gauge in which the fiducial universe 
has the same spherical symmetry as the LTB background. Next, in Section 4, we 
carefully examine how the polarisation evolves from the last scattering surface to 
the present time in the LTB spacetime with the help of the collisionless Boltzmann 
equation. In particular, we show that the temperature and polarisation anisotropy 
of CMB in the LTB model can be expressed by the same formula in terms of the 
shift vector as in the FLRW universe background. Putting these results together, 
in Section 5, we give general formulas for the CMB temperature and polarisation 
anisotropy produced by the void inhomogeneity and for the correlations among the 
harmonic components of the temperature, the E-mode polarisation and the B-mode 
polarisation. Finally, Section 6 is devoted to summary and discussions. 
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§2. CMB Polarisation Basics 

2.1. How to Represent Polarisations 

First of all, we explain the standard method to represent the polarisation of ra- 
diations. Let us consider a quasi-monochromatic plane electromagnetic wave prop- 
agating toward an observer, and take an orthonormal xy-basis that is orthogonal to 
the wave propagation direction. Then, the electric field of the wave is represented 
as E = ExBx + EyBy, with E^ = ax sin(a;t — ex) and Ey = ay sm{ujt — €y). 

In this setup, one can define parameters that represent polarisation as fol- 
lows: / := (a^) -I- {al), Q := (a^) - (a^), U := (2aj/aa; cos(€j; - Cx)), and V := 
{2ayax sm{ey — Cx))- These are called the Stokes parameters. Physically, I represents 
intensity (temperature), Q and U represent linear polarisation, and V represents 
circular polarisation. We ignore V because circular polarisation is never generated 
by Thomson scattering in the early universe. 

When the orthonormal basis is rotated in the wave plane, Q and U are linearly 
transformed. If we introduce the matrix P defined by 



-Q 

this transformation can be expressed as 

{e'x,e'y) = {ex,ey)R P' = R-^PR, (2-2) 

where R is a, two-dimensional rotation matrix, because the components of P can be 
written 

Pab = 4 {EaEh) - I6ab, (a, b = X, y). (2-3) 
Let us introduce the complex null basis defined by 

m:= ^{ex+iCy), ffi := ^{bx - iCy) . (2-4) 

Then, for the rotation of angle a represented by R{a), m transforms as 

m' = e^"m. (2-5) 

Hence, the complex quantity A by 

A:=Q + iU = m^m^Pah, (2-6) 

transforms as 

A' = e^'^^A. (2-7) 

Thus, A follows a simpler transformation law than Q and U ^ but still depends on 
the choice of the basis. This implies that we have to specify the basis in order to 
give a definite meaning to the values of Q and U. We will give such a specification 
below. 
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2.2. Polarisation Distribution Patterns 

Up to this point, we have considered a wave propagating only in one direction, 
whose polarisation can be described by the Stokes parameters {I,Q,U). In real 
observations, this set of parameters is measmed for photons of each direction, and 
the result is represented by three functions on the sky, /(0obs)i '3(^obs) ^^i^d ?7(^obs)) 
where 0obs represents the position on the sky. We define n := 0obs = (^obsj ¥'obs)- 
We often omit the subscript "obs" for brevity if no confusion occurs. 

Let denote the intensity (temperature) fluctuation around the sky average 
('2.725K'): 0{n) := AI{n)/{4I). This distribution is expanded by spherical har- 
monic functions as 

0{n) = J2&TYl^{n). (2-8) 

£m 

Inversely, 

0T = J dn0{n)Yp*{n) . (2-9) 

Next, let us turn to the polarisation distribution. In order to define the Stokes 
parameters Q{n) and U{n) for each direction n, we use the canonical orthonormal 
basis on the unit sphere with respect to the angular coordinates, {e0,e6i}. Then, 
as in the previous subsection, we can define the polarisation tensor P(n) on each 
direction n from Q and U, which can be expressed in terms of the complex Stokes 
parameters on the sphere, +2^(n) and -2A{n) = +2^(n), as 

Pabin) := +2A{n)rnamb + -2A{n)mamb , (2T0) 

where the subscripts a and b run over 9 and tp {6 = 1, (p = 2). These complex 
Stokes parameters transform under the angular coordinate transformation induced 
by a rotation of the sphere represented by a 3-matrix O as 

±2A'{n) = e±2^°o(")±2^(0-'n), (2-11) 

where aoin) is the rotation angle of the orthogonal frame induced by the trans- 
formation O. Functions on the unit sphere that transform in this way are called 
functions of spin-weight ±2. In general, these functions can be expanded in terms 
of spherical harmonics with spin- weight ±2 as 

±2^(n) = J2 ±2AT±2Yr{n) . (2-12) 
Here, the spin- weighted spherical harmonics are defined as^^^^'^^ 



^^" = yf^^"^^'^"'^"' -^^^^0; (2-14) 
sYp = , ^ < |s| ; (2-15) 



5 



where 



8, = -(sm«r(| + ^|){(sin«rM. (2-16) 
8, = -(si„«)- (A (2.17) 



or to be explicitly by 



{£ + m)l{i - my.{2i + 1) ^. 



+ (-l)^"'^"'e^'"^cot2'-+^-'"^ (2-18) 

Note that qY^ = Y^. These functions transform under the rotation O as 

,y/"(0-^n) = e-'"^°^''\Yf {n)Df'^{0), (2-19) 

where DJ^'^' {O) is the £-th irreducible representation of the 3-dimensional rotation 
group. From this it follows that the expansion coefficients ±2^™ transform under the 
3-rotation in the same way as the harmonic expansion coefficients for the temperature 
anisotropy 0J^: 

±2^7' = Df''^{0)±2Af. (2-20) 

Hence, we can define frame-independent harmonic expansion coefficients for the po- 
larisation distribution by 

Ef := ^{+2 AT + -2 AT) , BT := ^(+2^r - -2AT) • (2-21) 

These represent the gradient ( "E-mode" ) and rotational ( "B-mode" ) components of 
the polarisation field, respectively. 

2.3. Power Spectra from Primordial Fluctuations 

When the CMB temperature and polarisation anisotropics come only from pri- 
mordial fluctuations, without any secondary effects such as gravitational lensing 
effects, the power spectra are defined in terms of the correlation functions as 

{XiT*X2T ) = Cf^^^dee'Smm' , (2-22) 

or 

Cf'""' = E ' (2-23) 

m=—i 

with Xi, X2 e {0, E, B} and Xi, X2 G {T, E, B}. Here, the angle brackets 
represent an ensemble average over initial conditions; this average can be replaced 
in calculations with an average over space when the corresponding angular scale is 
much smaller than the observed region size. If physics and the initial condition are 
invariant under a parity inversion, we have Cj^ = Cf^ = 0. 
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§3. Gravitational Lensing Effects in the LTB model 

Inhomogeneous gravitational fields produce two effects on photon propagation. 
The first is a bending of its trajectory, and the second is the change of the pho- 
ton energy in addition to the standard redshift by cosmic expansion. The latter 
is the so-called Sachs- Wolfe effect, which we do not consider in this article. Intu- 
itively speaking, as far as CMB measurements by a fixed observer are concerned, 
the former — called 'shear field effect' — can be further divided into two parts: (i) the 
change of the photon direction in the sky and (ii) the displacement of the intersection 
sphere of the past light cone and the last scattering surface in the direction perpen- 
dicular to this sphere. In order to give a definite meaning to this distinction, we 
need to introduce some reference FLRW model to define 'unperturbed' photon tra- 
jectories and past light cones. However, this procedure introduces the gauge freedom 
corresponding to the mapping between the real universe and the reference model, 
and thus make that distinction obscure. In fact, for the FLRW model with small 
perturbations, the displacement of the last scattering sphere can be set to be zero by 
an appropriate gauge choice, and in this gauge, the shear field effect can be repre- 
sented only in terms of (i), namely, the 'gravitational lensing effect'. In the local void 
model, it is not so certain whether the same argument holds when the non-linearity 
of inhomogeneities is large. In the present article, we simply assume that the shift 
of the last scattering point in the direction normal to the last scattering sphere can 
be set to zero by a gauge choice. 

Under this assumption, the gravitational lensing effect on the CMB anisotropy 
can be simply determined by the two-dimensional shift vector 66ohs on the sky rep- 
resenting the difference between the observed direction of a photon and its initial 
direction on the last scattering sphere. 

3.1. Null Geodesies in the LTB Model 

Thus, the investigation of the gravitational lensing effect of a local void on CMB 

is reduced to determine 50obs as a function of the photon direction. For that, we 
have to solve the null geodesic equation in the LTB model, whose metric can be 
written 

ds'^ = -dt^ + S'^dx^ + r'^{de'^ + sin^ Od^p'^). (3-1) 

Here r is a function of t and x, and S is written in terms of r and the curvature 
function k{x) as 

S = r'/{\-k{x)x^fl\ (3-2) 

In terms of the photon four-momentum = dx^/d\ with affine parameter A, 
the geodesic equation can be written as dp^/dX = —Flfpp'^p^. Because of the spherical 
symmetry, this set of equations can be reduced to the coupled ODEs for w, and 
p± defined by 

pt = u;, p>' = nu;/S, pl=u;\l-fi^), (3-3) 



where 



P± ■.= r{{py + {p^f sin' ey/'. 



(3-4) 



(a) (b) 

Fig. 1. (a) Photon propagation in the LTB model. Each trajectory is contained in the unique 
two-plane passing through the center O, an off-center observer P, and the corresponding last 
scattering point Q. Without loss of generality, the two-plane with i/s = 0, vr is selected as that 
unique two-plane. The solid curve (red) represents the photon trajectory 7, while the dashed- 
dotted line (blue) represents the reference radial null geodesic 70. The solid circle represents 
the last scattering sphere for the off-center observer P, while the dashed circle represents that 
for the observer at the center O. (b) The shift vector SOohs- Here, only the 9 component SOohs 
is nonzero because of the above selection of the two-plane. 

Note that /U represents the cosine of the angle between the propagation direction of 
the photon and the radial direction of the photon position from the center. 

Without loss of generality, we can assume that the photon propagates on the 2- 
plane with (/? = 0, vr, and therefore p*^ = 0. Then, the geodesic equations are reduced 
to the set of four ODEs for uj{t), xi't)^ ^(t), and 

(3-5a) 
(3-5b) 

(3-5c) 

(3-5d) 

where ^ = (1 — k{x)x^y^'^ ■ Note that (3-5b) and (3-5d) form a closed set of ODEs 
for xit) and ii{t), and uj and 9 can be determined from each solution for and 
fi{t) by simply integrating knowing functions, using (3-5a) and (3-5c). 

We illustrate our setup in Figure la. Let O denote the center of the void and 
P the position of the observer in the comoving coordinate chart. We will take the 
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z-axis to run through O and P in order. Let us consider a radial null ray 70 with the 
angle ^obs relative to OP and a null ray 7 that passes through P with = — cos ^obs- 
Then, if D = OP is small, 7 stays close to 70 until the last scattering surface at 
t = tis. Hence, its behaviour can be determined by solving the linear perturbation 
equation obtained from the above ODEs for x (3-5b) and (3-5d). 

One subtle point of this approach is that /x cannot be treated perturbatively, 
because = — 1 for 70, but fi changes largely for 7 around P. This difficulty can be 
avoided by using the variables b and c defined by 

6 = X\/r^^, c = xiJ', (3-6) 

in stead of x and In fact, the geodesic equations can be written in terms of b and 
c as 

b = ^abc — I3b(? , (3-7a) 
c = i - ^ab^ + pb^c , (3-7b) 



S 



where a(t,x) and P(t,x) are 




" = ^(^-7) ' (^ = -2[--'^] ■ (3-8) 

Note that these are regular at x = in general and vanish for a spatially homogeneous 
background. Prom these, it follows 

The perturbation equation of the geodesic equation up to the first order in b 
now reads 

b=-Xb, (3- 10a) 

S' 

6c=-^Sc, (3- 10b) 

where X{t) is the function on the central past light cone x = Xoi't)^ 

^(t) = (xC« + x'/3)x=xoW • (3-11) 
These can be easily solved to yield 

b{t) = 6(io)e^(*'*°) , (3-12) 

where 

Y{a,b) = I dtX{t) . (3-13) 
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3.2. The Shift Vector J^obs 

In order to estimate the shift vector dOohs (Fig- lb), we need to calculate 6{t). 
We can assume that 6{to) = without loss of generality. In addition, assuming that 
d{t) monotonically decreases from ti to to, we choose the minus sign in (3-5c). Then, 
9{t) can be written 



dt x/r 



Jto XT Jto bl + ic/b)^ 
•*i {c/bydt x/r 



-I 



I, l + {c/byi/S-x^^a + x^cP 

(cmti) .tl ^c/bydt 



(c/&)(to) ^ + {c/by 



n {c/bydt f 
k ^ + {c/bf\ 



x/r 



+ [c/bf \ l/S - xHcc + X^cl5 



-'obs 



tan — 

X 



Jta I X^ r 



+ XP 



(3-14) 



where we have set c = — x that holds on the radial null ray. Further, the integral / 
in the last expression can be deformed with the help of a partial integration as 



Jto I X r 

Jto I V J Xdt 



1) 



^ bjti) - bjto) 

X{ti) 

Hence, we obtain 



+ 6(tc 



to 



(3-15) 



^(*ls) = ^obs - ^(*o) 

Xls 



*o dt / 1 _ X 



tls 



5 r 



(3-16) 



Thus, by eliminating the part that survives in the spatially homogeneous limit, we 
find that the shift in the angular direction of the null geodesic due to the void shear 
with respect to an observer P at distance D from the symmetry center O is given by 



59, 



obs 



-i^sin^obs r i^l - + ^ - ^e^^*'*")) , (3-17) 
4. X^ \r S r J 



where ^obs is the angle of the null geodesic direction with respect to the observer 
direction OP. This angular shift vector is used in Section 5. 
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§4. CMB Polarisation in the LTB model 

4.1. Flux Intensity Tensor 

In the Lorentz gauge, the free electromagnetic potential can be written in 
terms of the creation and annihilation operators as 

d^k 1 

p 

where ep^{k) is the polarisation basis satisfying 



/j3t, 1 



Note that the addition of vectors proportional to to Cp has no physical signifi- 
cance, because it corresponds to a gauge transformation and produces no physical 
effect in the exact quantum formulation. Up and satisfy the standard relativistic 
commutation relations 

[ap{k),a,{k')] = , [ap{k), a,(fc')^] = {27, f2u6pg6^k - k') . (4-3) 

The electric field E and the magnetic field B are 

Y,{ep{k)ap{k)e^'-^ + e;(fc)a,(fc)te-'=-) , (4.4a) 
(1^^^ X ^(ep(fc)a,(fc)e^^- + e;(fc)a,(fc)te-'^-) , (4-4b) 



B 

where 



kj 

^pj ~ (^pj ) (4'^) 

which satisfies 

fc • ep = , €*■ eq = Spq . (4-6) 

Now, let us define the measured components of the electric field, Sp, in terms of 
the sensitivity function W{x) and the detector polarisation basis €° as 

£p = J d^xW{x)e; ■ E{to,x) . (4-7) 

Then, for the free field, £p can be expressed as 
d^k 1 



" ' (27r)3 2 



aq(fe)(e° • eg{k))W{k)e-'''^ + c.c. 



(4-8) 



where 



W{k) = j d^xe'^ '^Wix) . (4-9) 
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Now, assume that 

(ap(fe)ag(fe')) = , (ap(fe)ta,(fe')) = 2{27rf ppg{k)d^k - k') . (4-10) 
Then, the observed correlation of the electric fields can be written 

, (4-lla) 

(: S,S, :) = I -^\W(k)\\-i)p[,j]{k) , (4-llb) 

where 

p,,(fc) = ^e;,(/e)e,,(fc)pp,(fc) , (4-12) 
p,g 

and £p is obtained from £p by advancing the phase by 7r/2 for each mode. Note 
that the Stokes parameters are 

p 

U = 2ere^Vfe) , V = (-2i)ere^Vfo-] • (4-13) 

Thus, Pij{k) provides a polarisation-basis- independent description of the radiation 
field polarisation and intensity. We call pij and ppq the flux density tensor and the 
fiux polarisation matrix, respectively. 

4.2. Polarised Boltzmann Equation 

In a curved spacetime, the above mode functions for the expansion of the elec- 
tromagnetic fields should be replaced by corresponding vector fields satisfying 

V^F^, = . (4-14) 

Accordingly, it is rather difficult to treat wavefunctions and polarisation vectors 
independently. However, such a treatment is allowed for modes for which the WKB 
approximation is good. For such modes, the mode function can be written 

A^{x) = a^(x)e^^(-) , (4-15) 

where for k := VS, a^{x) and S{x) satisfy 

k:=VS =^ k-k^O ^ Vkk ^ , (4-16a) 

Vkaf, = -^aSa^^O . (4-16b) 

Hence, by generalising the polarisation basis to spacetime-dependent vectors e'^{k,x) 
such that 

Vke'^ik, x) = 0, k^e^^ik, x) = , (4-17) 
the fiux density tensor pij can be generalised to 
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This tensor is independent of the polarisation basis and satisfies the generalised 
Boltzmann equation 

where C^^ is the collision term. 

4.3. The Initial Condition at the Last Scattering Surface 

On the last scattering surface, up to the linear order in perturbations, Jp^,^ can 
be Fourier decomposed into the contribution of each perturbation mode with the 
wave vector K as 

x,k) = j d^iife^^-pW {K- k) . (4-20) 
In this situation, it is customary to adopt the following polarisation basis: 

ep{k) = ^ep{k) ; (4-21) 

ei = = ik-{k- k)k) , (4-22a) 

62 = =k X k , (4-22b) 

^l-{k.k)^ 

where a is the scale factor*^ and 

In the spherical coordinates corresponding to the Cartesian coordinates in which 

K = (0, 0, 1) , k = (sin 9 cos ip, sin 9 sin ip, cos 9) , (4-24) 

the polarisation basis has the following components: 

€i = (cos 9 cos if, cos 9 sin ip, — sin ^) , (4-25a) 
€2 = (— sin f/:;, cos 0) . (4-25b) 

For this choice of the polarisation basis, because all vector-like quantities for a per- 
turbation is parallel to K, p^g defined by 

pS^=er<pg(^;fc) (4-26) 

depends only on oj and n = cos 9 = K ■ k: 

p«=pjy(i^;u;,/x). (4-27) 

Based on this setup, wc can solve the Boltzmann equation (4-19). This can be 
done following the standard procedure used in the case of the FLRW universe^^^^^^ 
and finally we obtain the initial condition at the last scattering surface. 

*^ The early universe can be treated as the FLRW universe plus a perturbation as described 
below. 
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(a) (b) 

Fig. 2. (a) Propagation of the polarisation basis in a general perturbed LTB spacetime. The polari- 
sation basis eq{tis, x{xo, ko), k{xo, feo)) at Q(fco) on the last scattering surface t = tis is parallelly 
propagated along the null geodesic 7(P, feo) to P on the present-time hypersurface t = to and 
named , which should be compared to the polarisation basis ep{to,xo, feo) at P. (b) The global 
polarisation basis in the unperturbed LTB background. The polarisation basis ep{z,k) at the 
present void center O, together with fe, is transported to an arbitrary point P on the hypersur- 
face t = to along the radial line J7 — const, parallelly with respect to gij{to,x) and becomes 
ep(P,fe). Meanwhile, e'p{z,k) at O is transported back to the symmetry center O' at t = tis 
along the timelike path corresponding to the symmetry center to define the basis at O', which is 
consecutively extended to an arbitrary point on the last scattering surface t — tia by the parallel 
transport with respect to gij{tis, x) and becomes ep(tis, a;, fe). 

4.4. Propagation After Last Scattering 

Let us work in the synchronous gauge in which 

ds"^ = -dt^ + gij{t, x)dx'dx^ , (4-28) 

where cnj approaches a spatially homogeneous and isotropic metric gij in the early 
universe: 

gij = 9ij + Sgij . (4-29) 

We take the constant-time surfaces so that the last scattering surface is represented 
hy t = tig. Note that the LTB spacetime g'^^^ also belongs to this class: 

9ij = 9ij + SuTBgij ■ (4-30) 

In most part, we take the spatially flat FLRW solution as the background g^j^. 

Now, we consider the measurement of the CMB polarisation at the space- 
time point P = (to,iCo). The null geodesic 7(P,fco) passing through P with the 
four-momentum = (uj{kQ),kQ) intersects with the last scattering surface at the 
point Q(/co) with the space coordinates x = x{xo, ko) and the four-momentum 
k = k{xo, ko). Let us represent polarisation bases at each point on the hypersurfaces 
t = to and t = tis as ep{to, x, k) and ep(tis, x, k), respectively. Then, the flux polari- 
sation matrix at P is expressed in terms of the corresponding quantities on the last 
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scattering surface as 

p{F,ko) = Cp{tis,x,k)& , (4-31) 

where x = x{xo,ko) and k = k{xo,ko) are understood. The matrix C is defined as 
follows. First, we parallelly propagate the polarisation basis eq{t\s, x{xo, ko), k{xo, ko)) 
at Q(A;o) on the last scattering surface along the null geodesic 7(P, ko) to P (Fig. 2a). 
Let Cq denote this basis at P. Then, 

Cpq = ep{to,xo,ko)-ef . (4-32) 

When the universe is well-described by the FLRW model at and before the last 
scattering, p{t\s,x,k) is approximately isotropic: 

ppqitls, X, k) = ^Spql (^-^^ + dppqitis, X, k) . (4-33) 

Then, from the identity 

J2CpqC;, = Sqr, (4-34) 
P 

we obtain 

Ppq(P, ko) = ^Spql + CprC*^dprs{t\s, X, k) . (4-35) 

The first term on the right-hand side of this equation is only relevant to the temper- 
ature and the polarisation comes only from the second term. In this formulation, it 
is understood that 5prs{tis,x,k) is for a polarisation basis that does not depend on 
the mode wave number. 

Now, we show that Cpq is close to the unit matrix in an appropriate global 
polarisation basis(Fig. 2b). First, we define the polarisation basis at the center O of 
the void at present as 

ei(z, k) = ^ (z-{k- z)k) , (4-36a) 

^Jl-{z.k)^^ 

62(2, k) = ^ --Z X k , (4-36b) 

^^\-{z.kY 

where z is a unit spacelike vector and k = k/\k\ as before. Later, z is taken to 
be the direction of the observer from O. Next, we transport this pair and k'^ to an 
arbitrary point P on the hypersurface t = to along each radial line parallelly with 
respect to the space metric gij{to,x) to define the polarisation basis at P with respect 
to k, ep^(P,k). In this way, we can define a polarisation basis everywhere on the 
hypersurface t = to. 

Next, we parallelly transport the basis at O back to the symmetry center O' 
at t = tis along the timelike path corresponding to the symmetry center to define 
the basis at O'. Then, extend it to an arbitrary point on the last scattering surface 
t = tis as on the hypersurface t = to. When the spacetime is exactly spherically 
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symmetric, we can easily relate this basis to the mode-dependent basis introduced 
in the previous section, if we neglect the small spatial inhomogeneity of the LTB 
universe at t = tig. Because the matrix Cpq that we are calculating is multiplied 
by a perturbation 6ppq in (4-35), and wc know 5ppq only in the first order w.r.t. 
perturbations at last scattering, we can neglect the small inhomogeneity of the LTB 
universe at the last scattering in calculating Cpq. 

Note that when the polarisation basis €p{k) is given for a background space 
metric gij, we can uniquely determine the corresponding basis ep{k) for gij = gij+Sgij 
by the requirement 

ep-eq = Spq , €p ■ fc = (4-37) 
as ^ 

^^Pj = -^hjeepmSg^"" , (4-38) 

where 

hij = 9ij - kikj . (4-39) 

Hence, in order to determine the matrix Cpq, wc only have to calculate the propa- 
gation of the polarisation basis for the background LTB universe neglecting small 
inhomogeneities. 

Now, we show that the transfer matrix Cpq for the polarisation basis can be 
well approximated by the unit matrix for the global polarisation basis introduced 
above. Because we are only interested in the shear effect of the LTB geometry, we 
can calculate this matrix by parallelly propagating the unperturbed basis e^^ in the 
background LTB universe neglecting additional perturbations. Now, let us consider 
a null geodesic 70 passing through O at t = Iq and a null geodesic 7 passing through 
an off-center observer at P close to O. Then, it is easy to confirm that for the null 
vector k parallel to 70, the corresponding polarisation basis e!p^{t,k) at t = to and 
t = t\s defined above are parallelly related. 

Let ep^{t,z) be a polarisation basis parallelly propagated along a family of null 
geodesies with the parameter z close to 70 among which 7 is contained. Then, along 
7, we obtain 

Vfe[e; • (Vze;')] = e'q ■ R{k, Z)e'p , (4-40) 

where V z refers to the covariant derivative along the deviation vector Z. This can 
be written 

|« • Vfee;') = e'q'e'p'Rje^m Sx^ , (4-41) 

where 5x is the short distance along Z. Here, the polarisation basis e'p on 7 can be 
written explicitly in the spherical coordinates for the LTB metric as 

e'^ = --de, e'^ = -^d^. (4-42) 

r r sm ti 

Further, the non-vanishing components of the Riemann tensor for the LTB metric 
are 

Rt^tx = -SS , RtAt"" = ~5a. RtAx"" = ^ 5 A , (4-43a) 
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XX 


= SS , 


pt 

AB 


= rrjAB , 






pX 


S 

~ S ' 


pA _ 








px 

XX 


S' 

~ S ' 


pX 
^ AB 


rr' 

= --g2^AB , 


T^A cA 
XB - —^B > 





i?,^/ = ^ - - (^-J , (4.43b) 

Rab^'' = ^ I'^' + 1 - (^)'} (^aSE - S^S^b) . (4-43C) 

From these, it follows that the right-hand side of (4-41) vanishes. Further, from the 
structure of the non-vanishing components of the Christoffel symbol for the LTB 
metric, 

(4-44a) 

(4-44b) 

(4-44c) 
{jab is the S'^ metric) 

we find 

ViW = DiV^ + KjV^ , (4-45) 

where y is a four- vector, Di is the spatial covariant derivative with respect to gij 
with constant t and K^^ is the extrinsic curvature. Since the temporal component 
vanishes for a vector along constant-time hypersurfaces and the spatial covariant 
derivative with respect to gij along Sx vanishes for a vector parallelly transported 
with respect to gij along Sx, this implies that 

4 • Vfec; = (4-46) 

for t = to and t = t\s. Therefore, we have found that we can set C to be the unit 

matrix for our choice of the polarisation basis, and that the observed polarisation 
tensor can be expressed in terms of the perturbation of the polarisation tensor at 
the last scattering simply as 

Ppq{t0, Xo, ko) = ^Spq (^-^^ + 5ppq{tls, X, fc) . (4-47) 

Thus, the expressions for the temperature and polarisation can be calculated as in 
the FLRW-universe case. 

Finally, we note that we have to change the polarisation basis from the global 
basis to a mode-dependent one in order to express dppq{t\s, x,k) in the above formula 
in terms of the mode-dependent expression for the initial condition, p(^) in §4.3. 

From the spatial flatness at last scattering and the spherical symmetry of the 
background, it follows that the global polarisation basis on the last scattering surface 
defined above can be written as 

4 = , (4-48) 

where e^'s are vectors that have the same expression in the coordinate system in 
which the FLRW background metric at the last scattering is expressed as 

ds'^ = -df + a{tfdx ■ dx . (4-49) 
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It is easy to find the relation between this and ep{K, k). It is given by the rotation 
matrix 

Apg := e'p{tis, X, k) ■ eq{K, k) = RpqiOi^) (4-50) 
with the angle 9\s{K,k) satisfying 

K-z-(k-k)(i-k) . , K-(kxz) 

cos 6'is = — , sin 6'is — 



'1-{K- kf^Jl -{z-kf ^^-{K- kf^Jl - (i • kf 

(4-51) 

The initial condition for 5ppq at the last scattering surface is thus given by 

5ppq{t = tis, x,k) = J d^Ke'^ ^'AiK, k)p^'^'^ (K, u, ii)A{K, k)^ . (4-52) 

§5. Results 

5.1. Temperature 

We denote the unlensed and lensed temperature anisotropics by and 0', re- 
spectively. By the shift vector SO, we define the direction n'(n) := n + 50, where 
n := —ko is the direction to which the observer looks. Prom (3-17), the components 
of the shift vector are 69 = DF sinO and Stp = 0, where 

Therefore, if D = 0, then n' = n. With these definitions, we have 

0'(n) = 0{n') . (5-2) 
The l.h.s. of (5-2) can be expanded as 

0'{n) = Y,O7Yr'{n), (5-3) 

im 

while the Taylor expansion of the r.h.s. gives 

0{n') = 0{n + 50) = 0{n) + 50 ■ V0{n) + 0{D'^) . (5-4) 
Hence, 50J^ := 0'^ — 0^ can be calculated up to the first order in I? as 

5&f = Dr j dnY['*{n){p'^ - 1)9^ ^ 0fY^'{n) 

e',m' 



Dr 



4^2 _ 1 «-i ^ ' '\j A{1 + 1)2-1 
where we have used the identity^''') 



(5-5) 



.2 d 



ip' - l)—Yt = iCT+^YHl, + l)CryA , (5-6) 
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with 



5.2. Polarisation 



As in the case of the temperature anisotropy, the lensed polarisation tensor 
P^j(n) and the unlensed polarisation tensor Pabin) are related by 

PM = Pab{n'). (5-8) 

Accordingly, we obtain the relation 

E'T±iB'T=J2iEf±iBF') I dn±2Yr{n){l^'-l)d,±2Yjf'in) . (5-9) 



',m' 



Using the formula^''-* 

(/x^ - i)§^sYr = isCT^^sYr^^ - + \)sCTsYr-^ + j^y^r , (s-io) 

with 



" V ^2(4^2 _ 1) ' (5-11) 



we obtain 



.r,m _ ((^ + l)2-m2)((£ + l)2-4) 



(£2-m2)(^2_4) 



((^ + l)2-m2)((£ + l)2-4) 

--^^W (^ + i)2(4(£ + i)2_i) (^ + 2)5m 



(^2_^2)(^2_4) 

Prom these, we sec that even if there exists no B-modc in the unlensed polarisation 
anisotropy, an off-center observer with Z) 7^ detects non- vanishing B-modes in the 
lensed polarisation anisotropy. 



5.3. Correlations 

Prom these relations and the assumptions on the correlations for the unlensed 

-fTB _ ^e: 

'I - 



initial anisotropics, (2-22) and Cj^ = Cf^ = 0, we immediately obtain the following 
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formulas for the correlations among lensed anisotropies of the CMB temperature and 
polarisation: 



+ 1)2 — m? 



+ + 1)2 _ 1 {^g/ - + 'i)Cl+\}5t,e-iK 



£2 



2 



+ Dr\l ^-^^^{-{l + l)Cr ^ _ i)qTT}5^^^,^^5^^, (5.13a) 



(£+l)2_^2 (£ + l)2-4 ^^TE r^^9lrTElx . 



(5-13b) 

, n r / ((^+l)^-m2)((£ + l)2-4) EE , .^^EE . . . 

+ (4(£+l)2-l)(^ + l)2 i^C-, - + 2)Q+i} V-K^W 



(5-13c) 

(©7*57') = -iDrj^^Cr5,,6mm' (5-13d) 

(£;7*i?7'> = -^^^^(^^(^1''' + CrW^ra' (5-13e) 
{B'^*B'^ ) = Cf^6u'Sjmn' (5-13f) 

These expressions show that an off-center observer detects non- vanishing T — B 
and E — B correlations. This is what was expected, but these correlations have one 
non-trivial significant feature; they are nonzero only for the same ^, while [E'f'E'"} ) 
are nonzero for = £ ± 1. This feature can be used as a decisive signal showing the 
existence of a local void in future B-mode measurement experiments with very high 
precision. 

5.4. Estimation 

Although the main purpose of the present paper is to derive general formulas 
for the temperature and polarisation anisotropies in the LTB model, we give some 
numerical estimations for F in order to get a quantitative idea on the lensing effect 
of a local void. 

We consider two models characterized by the curvature parameter function 
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Table I. Numerical estimates of F for the AA model and the modified AA model. 



Void radius 
-ffoxo 


Wall width 
HoAx 


Value of r/Ho 


AA model 


Modified AA model 


0.0235 


0.020 


-3.62 


1.68 


0.235 


0.20 


-0.09 


0.30 


0.470 


0.40 


0.125 


0.303 



Density contrast parameter a = 0.9 



i^kix) the density parameter function f2m{x) defined by 

r2,(x) = --^fc(x), ^m{x) = ^-^^, (5-14) 

where k{x) is the curvature function appearing in (3-2), and M(x) represents the 

mass inside the sphere of the coordinate radius x a-t present t = to- 

The first model (AA model) is the void model studied by Alnes, Amarzguioui 
and Gr0n in 11) and 12). This model is defined by 

^kix) = f (l - tanh , (5-15a) 

= 1 - ^kix)- (5-15b) 

The second condition fixes the gauge freedom in the choice of x- This model ap- 
proaches the Einstein-de Sitter model at infinity outside the void, and the ratio of 
the expansion rates of the universe at present at infinity and at the center, H^/Hq 
is given by 2/{SHQto), where to is the age of the universe at the center. 

The second model (modified AA model) is a modification of the AA model and 
defined by 

^kix) = f (l - tanh ^^^p^) , (5-16a) 
Hmix) = 1 - ^fc(x)- (5-16b) 

The integrand for F in (5-1) depends on the second x-derivative of metric coefficients. 
Because this model is smooth at the center, the integrand is finite at the center. In 
contrast, the AA model has a cusp singularity in the density and curvature. Hence, 
the corresponding integrand has a kind of (5-function type singularity. Although we 
can obtain a finite value for F by neglecting this contribution, it comes up when we 
replace the model by a smoothed one. This is the reason why we considered this 
modified AA model. 

In Table I, we give the results of numerical estimations for three model parame- 
ters for each model. Roughly speaking, F/Hq is of order unity, hence 69 = 0{HqD). 
However, the exact magnitude and sign are quite sensitive to the void shape and size 
parameters. This should be contrasted with the dipole anisotropy, which is rather 
insensitive to the void size and shape. More detailed analysis is under investigation. 
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§6. Summary and Discussion 

In this paper, we have developed a formulation to calculate the gravitational 
lensing effects on the CMB temperature and polarisation for an observer close to the 
center of a spherically symmetric void described by the LTB model. In particular, we 
have derived explicit expressions for the correlations among the anisotropies of tem- 
perature and polarisation induced by gravitational lensing in terms of an integration 
of known geometrical quantities along the central past light cone. 

With the helps of these formulas, we have found that for an off-center observer 
in the local void, there appear nonzero correlations between T and B and between 
E and B that are diagonal in the harmonic coefficient expression in the leading 
order with respect to the observer offset distance. Similar correlations arise if there 
exists a quintessence-type axionic field with mass in the range from 10~^^eV to 
10~^^eV, but in this case the magnitudes of correlations have different dependence 
on i. Hence, if the correlations suggested by our result are detected in future B- 
mode measurement experiments with high precision, they would provide a clear 
signal showing the existence of a local void. 

We have also given some preliminary numerical estimations of the gravitational 
lensing effect. The results indicate that the B-mode amplitude produced by the 
lensing is around 10~^ times that of E-modes if we take into account the constraint 
on the off-center distance on the observer from the dipole anisotropy of the CMB 
temperature. -"^^^ This is because the effect is proportional to the distance of the ob- 
server to the symmetry center. Thus, we need next generation B-mode experiments 
to use this effect to test the viability of the local void model. However, the results 
also show that the gravitational lensing effect is very sensitive to the void profile 
and sizes. Hence, when it is detected, it is useful to specify a model. Because the 
lensing effect becomes larger for smaller voids in general, it may be also used to 
detect anomalously large voids, which are much smaller than the standard void size 
(~ IGpc) in the local void model but still statistically rare in the ylCDM model. 

In the present paper, we have only considered the local void model, but the 
method developed here can be applied to more general inhomogeneous models, such 
as realistic universe models with voids outside us. Such extensions are under inves- 
tigation. 
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